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We explore the extension of quantum cosmology outside the homogeneous approximation, 
using the formalism of loop quantum gravity. We introduce a model where some of the 
inhomogeneous degrees of freedom are present, providing a tool for describing general fluc- 
tuations of quantum geometry near the initial singularity. We show that the dynamical 
structure of the model reduces to that of loop quantum cosmology in the Born-Oppenheimer 
approximation. This result corroborates the assumptions that ground loop cosmology, sheds 
light on the physical and mathematical relation between loop cosmology and full loop quan- 
tum gravity, and on the nature of the cosmological approximation. Finally, we show that the 
non-graph-changing Hamiltonian constraint considered in the context of algebraic quantum 
gravity provides a viable effective dynamics within this approximation. 

I. INTRODUCTION 

Loop quantum cosmology (LQC) provides the most successful physical application of loop grav- 
ity, and one of the most promising avenues towards the possibility of an empirical test of quantum 
gravity [U [2]. (For a recent review of quantum cosmology, see [3].) In particular, the possibility 
of a fully consistent quantum description of "big-bang physics" , and the robustness of the bounce 
prediction, represent a clear advance in our understanding of quantum gravitational physics within 
this theoretical framework. This impressive success opens a number of physical and mathematical 
questions: (i) Can we include inhomogeneities? Inhomogeneities and their quantum fluctuations 
play a fundamental role in the currently fashionable cosmological scenario. Inhomogeneities can in 
principle be re-inserted at a later cosmological epoch, restricting the analysis of the Planck epoch 
to the sole homogeneous degrees of freedom — but is this approximation viable? After all, what 
is very interesting is precisely to understand the configuration of the full fluctuating quantum ge- 
ometry near the singularity itself, about which very little is known (see [1]). In other words: (ii) 
Can we describe the actual quantum state of the geometry near the initial singularity, beyond the 
homogeneous approximation? Perhaps this state could even teach us something directly about the 
emergence of the physical inhomogeneities of our universe, (iii) What is the true relation between 
full loop quantum gravity [5j |6] and LQC? The question has been addressed repeatedly p 7 ] and 
concerns have been raised on whether the two theories are truly consistent — some simple minded 
ways of interpreting their relation have even been recently rigorously proven incorrect [B] • 

We address here all these questions. We do so by first analyzing the nature of the approximation 
on which cosmology itself -classical or quantum- is based. This is neither a low-energy nor a high- 
energy approximation, since cosmology appears to describe well very large distance features of our 
universe as well as its behavior in much higher energy-density regimes. The analysis leads us to 
the idea that the full theory may be consistently expanded by adding degrees of freedom one by 
one, starting from the cosmological ones. Accordingly, we define an approximated dynamics of the 
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universe, inhomogeneous but truncated at a finite number of degrees of freedom, and we discuss its 
regime of validity. (For previous works on inhomogeneities in LQC, see [§1 \1U\ ITT]. See also |12j.) 
This approximation includes and extends conventional cosmology, without however including the 
full infinite dimensional field theory. We work in the compact case, which is conceptually simpler 
— that is, we assume that the topology of the spacial universe is that of a three-sphere (a possibility 
which is still compatible with observations, and, according to some [13], it is even favored by them). 
The approximation we take can be intuitively interpreted as a truncation of all degrees of freedom 
to a finite order in a multipolar expansion of the fields on the topological three-sphere. These 
degrees of freedom can be described using a fixed 3d compact triangulation A n , formed by n 
tetrahedra. 

The quantum kinematics of this system turns out to be described by the truncation of loop 
quantum gravity obtained by restricting the spin-network states to those based on a graph equal 
(or contained in) to the graph defined by the dual of A n . Within the approximation considered, 
the quantum dynamics can be described by the non-graph-changing version of the Hamiltonian 
constraint [J3] that has been recently considered in the context of algebraic quantum gravity JE] . 
Thus, non-graph-changing Hamiltonian constraint plays here the role of an effective dynamics, as 
originally suggested by Thiemann 1 . In this way, one can try to define a quantum cosmological 
model for any given triangulation A n , with a number of degrees of freedom that increases with the 
complexity of A n . 

We study here in detail the simplest nontrivial case, based on a triangulation A2 formed by the 
minimal triangulation of a 3-sphere: two tetrahedra glued along each face. The corresponding dual 
graph is formed by two nodes joined by four links. We show that the model is well defined and 
in particular the constraint algebra closes. We write the state space, the quantum operators and 
the Hamiltonian constraints of this model explicitly. The model represents a non homogeneous 
quantum universe, where, say, we do not consider just the overall spacial average of a scalar field, 
but also its dipole moment. This provides a well-defined "first step out of homogeneity" in loop 
quantum cosmology. We call this model a "dipole cosmology" , but the name should not be taken 
literally, as the gravitational degrees of freedom that are described are more "quadrupolar" that 
dipolar: the spatial universe is split into two hemispheres (represented by the two tetrahedra), 
separated by a closed surface E; but S splits in turn into four large surfaces (the four triangles 
bounding the tetrahedra) , whose areas provide degrees of freedom that can be roughly thought as 
capturing the geometry of a 3d ellipsoid. 

We then ask in which sense LQC is contained in the larger model. We argue that a proper way of 
addressing the problem is to interpret LQC as a first-order Born-Oppenheimer approximation. This 
is the approximation in which the effect of the inhomogeneities on the dynamics of the scale factor is 
small, as the effect of the electrons on the dynamics of the nuclei is small, in the Born-Oppenheimer 
approach to molecules [TE] • (For other utilizations of the Born-Oppenheimer approximation in this 
context see [E]-) We show concretely that taking the order-zero Born-Oppenheimer approximation 
of the A2 model yields precisely the structure of the LQC dynamics. In particular, we recover the 
characteristic structure of the LQC Wheeler-DeWitt equation, defined as a 3-terms finite-difference 
equation in the scale factor. This derivation provides a prototype for understanding how LQG is 
contained in full loop quantum gravity. In particular, we derive here the quantization of the LQC 
\i parameter directly from LQG, without need of an explicit recourse to the area gap argument PQ. 

The models defined here, and in particular the A2 dipole cosmology, are finite-dimensional 
quantum theories, that can be used to describe the inhomogeneous quantum geometry near the 
initial singularity and its quantum fluctuations, to any given arbitrary order. We leave the analysis 
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of these models to further investigations. 

The paper is organized as follows. In Section 2, we discuss the nature of the approximation we 
take and its viability. In Section 3, we define the quantum models for arbitrary A n . The A2 model 
is described in some detail in Section 4. The Born-Oppenheimer approximation is introduced in 
Section 5, where we show how the structure of LQC can be recovered. In Section 6 we conclude 
pointing out some implications of the results presented, in particular their relation with Regge 
calculus cosmology [18], and with the approximation used in computing background independent 
ra-point functions [19]. 

II. APPROXIMATIONS IN COSMOLOGY 

Modern cosmology was born with Einstein's 1917 seminal paper [20], where Einstein states 
the cosmological principle, according to which the dynamics of a homogeneous and isotropic space 
approximates the one of our real universe. A certain vagueness lingers around the precise status 
of this principle, presented under a variety of different lights in the literature. The question of 
whether the universe approaches homogeneity at large scale is an empirical question: the later 
evolution of observational cosmology appears to be corroborating it, as well as finding some pre- 
cise quantitative limits to large scale homogeneity. Einstein's principle is therefore a hypothesis, 
in the healthy tradition of hypothetical-deductive science. But what precisely is the hypothesis, 
and why is it relevant? The universe is obviously not homogeneous at every scale and not exactly 
homogeneous at any scale, except at its very largest scale, where it is homogeneous by definition. 
Why can we neglect the effect of the inhonomogeneities on the dynamics of this largest scale, in a 
nonlinear theory like general relativity, where all scales are coupled? In fact, this is precisely the 
hypothesis put forward by Einstein in 1917: that the universe happen to be in a state where the 
effect of the inhomogeneities on the dynamics of its largest scale, described by the scale factor, 
can be neglected in a first approximation. More generally, that the universe is in a configuration 
where the effect on the longest-wavelength of the interaction with the shorter- wavelengths is neg- 
ligible. The significance of this hypothesis is becoming particularly clear today, since a number of 
contrary hypotheses are being explored, such as, in particular, the very intriguing possibility that 
the measured cosmological constant could be -in full or in part- the result of the effect of these 
shorter wavelengths on the largest scale (see for instance [21]). In other words, the cosmological 
principle is the hypothesis that a certain approximation scheme is viable in general relativity, and 
that the universe happens to be in the regime where this approximation scheme is effective. 

What precisely is this approximation? One is tempted to say that it is simply a long-distance 
one: it is defined by cutting-off the modes with wavelength shorter than a certain size L. But this 
is imprecise, because L varies with the size itself of the scale factor, and can be also very small. 
The situation is easier to analyze in the context of a spatially closed universe; let us therefore 
assume in the following that we are in this context. Let a 3 (i) be the volume of the universe at 
the cosmological time t. Then the approximation on which cosmology is based is to neglect the 
dynamics of the wavelengths A shorter than a(t). But a{t) itself can be large as well as small. If 
n ~ j, the cosmological approximation is the first order term, n = 1, in an expansion in small n, 
which does not necessarily mean large A. 

Consider the next terms in the expansion, say n = 2,3,... and so on. An immediate conse- 
quence of the cosmological principle (or a natural extension of the same) is the hypothesis that 
the dynamics of the long wavelengths modes of the universe is only weakly affected by its highest 
modes. We can implement this expansion by approximating the geometry of the universe not just 
by a maximally symmetric space, as in standard cosmology, but rather by a geometry described 
by a finite number of degrees of freedom. These can be labeled by the elements of a triangulation 
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A n formed by a finite number n of tetrahedra. If a maximally symmetric space represents the 
gravitational degrees of freedom of the universe averaged over the largest possible scale, the de- 
grees of freedom of a fixed triangulation can be interpreted as a representation of the gravitational 
degrees of freedom of the universe, averaged over the large scales, up to a certain degree in a mode 
expansion. In the following, we consider the classical and quantum description of the dynamics of 
a model of the universe defined in this manner. That is, we construct an effective theory where 
wave lengths A < Aq ~ ^ are neglected. 



III. THE MODEL 
A. Classical theory 

Fix an (oriented) triangulation A n of a (topological) three-sphere, formed by n tetrahedra t 
glued by their triangles. We label the triangles with an index / ("/" for face) that runs from 1 to 
2n (the number if faces is twice the number of tetrahedra). Associate a group element Uf € SU{2) 
and a su(2) algebra element Ef to each oriented triangle /. We take the convention that to the 
face f~ l obtained inverting the orientation of / is associated the group element 

Uf-i = UJ 1 (1) 

and the algebra element 

E f -i = - Uj'EfUf. (2) 

If Tj, i = 1,2,3 is a basis in su(2), we write Ef = E l ^T{. We take Uf and Ef as phase space 
variables of a dynamical system, with the conventional Poisson brackets structure of a canonical 
lattice SU(2) Yang-Mills theory, that is 

{U f ,U f ,} = 0, (3) 
{E),U f/ } = 5 f f,T*Uf, (4) 
{E),E fl } = -bff.e^E). (5) 

In other words, the phase space is the cotangent bundle of SU(2) 2n with its natural symplec- 
tic structure. Let the dynamics of the system be defined by two sets of constraints: the gauge 
constraints 

Gt = Y,E f ~Q, (6) 

fet 

where the sum is over the four faces of the tetrahedron; and the Hamiltonian constraints 

C t = V t ~ l Yl T±[U fr E f ,E f ,]~0 (7) 
ff'et 

where the sum is over the couples of distinct faces at each tetrahedron, Uffi = UfUf 1 Uf 2 ...Uj, 1 
where Iff = {/, /i, fi, is the link of the oriented faces around the edge where / and /' 

join, and 
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V t=\ E T±[EfE f ,E r ]. (8) 
fff"& 

where the sum is over the four unordered triplets of distinct faces at the tetrahedron and {/, /', /"} 
has positive orientation. Notice that = Tr[EfEf/Ef»] because of ([6]). This concludes the 
definition of the dynamical systems we want to consider. In Section IV we study one of these 
systems in more detail, and we make sure the constraint algebra closes. We leave the analysis of 
the constraint algebra in the general case for future developments. 

This dynamical system can be interpreted as a cosmological approximation to the dynamics of 
the geometry of a closed universe. To see this, consider real Ashtekar fields A l a (x) and E ia {x), 
with their standard Poisson algebra (see for instance [5]), on a 3d surface £ with the S3 topology. 
(The index a is a 3d (abstract) tangent index.) Let A n be a triangulation of £ and A* a dual 
of the triangulation. We interpret Ut as the parallel transport of the Ashtekar connection A a {x) 
along the link e/ of A* dual to the triangle /, and Et as the flux $j of the Ashtekar's electric 
field E a {x) across the triangle /, parallel transported to the center of the tetrahedron. That is, 
Ef = Uf^fUJ 1 and Ef-i = —UJ^<& fU '/ 2 where Uf = UfUf 2 and Uf and Uf are the holonomies 
of the two segments ef , e/ 2 into which the face / cuts the link e / . Then the Poisson brackets of 
Uf and Ef defined in these manner turn out to be precisely ( 3|4|5 ). In particular, notice that the 



origin of Q5J) is the fact that is parallel transported to the center of each tetrahedron. (Notice 
also that ([5]) follows from (pjj), Q and the Jacobi identity.) 

The gauge constraint d6J) generates the correct internal gauge transformations on these vari- 
ables. 2 If the triangulation is sufficiently fi ne, ^ approximate the Ashtekar's (euclidean part of 
the) Hamiltonian constraint Tr[F a f,E a E b ] / V detE ~ 0, where F a i, is the curvature of A a . Notice the 
absence of the second term of the usual discretization Tr[(Uff — UZ,)EfEf] of the Hamiltonian 
constraint. The Ujj} is usually subtracted in order to subtract the first term in the small-curvature 
expansion U ~ exp J A ~ 1 + |a| 2 .F a b + O(|o| 4 >l 2 ); but the subtraction is not needed because this 
term does not contribute to Ct thanks to ^ . We do not know if this observation has already been 
made in the literature. An explicit calculation shows that Ct is real. 

Finally, Vt is (proportional to) the volume of the tetrahedron t and we call V = Vt the total 
volume of space. We take for convenience 8/37rG Newton-, the speed of light and h to be unit, and 
we choose the Immirzi parameter 7 = 1 for simplicity; what follows needs to be extended to the 
more interesting case of real 7 and full Hamiltonian constraint. 

The constraint ([T]) corresponds to the non- graph-changing version of the Hamiltonian constraint, 
as the one utilized in algebraic quantum gravity [15]. This is the only viable alternative in the 
present context, where we have reduced the degrees of freedom of the gravitational field to a fixed 
number. This version of the constraint approximates the classical Hamiltonian constraint if Uff, 
namely the parallel transport of A a along the loop a dual to the link Iff, approximates 11 + |a| 2 F a fc. 
It is important to notice that this happens not only if the length of the loop is small, but also for 
large loops if A a is small. Hence near flat spacetime the approximation can be good even for coarse 
triangulations. Misunderstanding of this fact has generated the erroneous idea that low-curvature 
spacetime needs to be approximated by fine triangulations. 

Alternatively, one can interpret the data (A n , Uf, Ef) as a description a piecewise flat Regge ge- 
ometry, where the curvature of the connection is concentrated on the edges of A. The quantity Uff 
gives then the curvature at the corresponding edge. However, observe this interpretation is slightly 
misleading here, since the variables of the model are better understood as macroscopic quantities 



2 In particular, writing G[X] := 2^2 t Tr[XtGt] where At G sit(2), the infinitesimal gauge transformation of Uf is 
SUf — {Uf,G[X]} = Xt t Uf — UfXt 2 where (ti,ta) are the two tetrahedra separated by /. 
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averaging over local degrees of freedom. Therefore the flatness of the individual tetrahedra does 
not need to be take literally. 

It's simple to couple a family of multifingered "clock" variables, one per node. The simplest 
choice [I] is an (ultra-local) scalar field 3 with a value cftt and conjugate momentum p^ t at each 
node, with the overall Hamiltonian constraint given by 

77 £ T *\- U ff' E f' E f] + W P *t~ °' (9) 
t fff& t 

where k, proportional to the Newton constant G, determines the matter-gravity coupling. The 
role of this field is double. First, it keeps track of evolution in a background-independent manner, 
namely it models a physical clock. Second, it represents in a simplified manner the matter content 
of the universe. Replacing this field with a more realistic description is viable here: ultralocality 
can be eliminated adding a difference term; while Yang-Mills and fermion fields have a particularly 
straightforward description in this language [5 J . 



B. Quantum theory 

The quantization of the model is immediate. Following what is done in lattice QCD, a quan- 
tum representation of the observable algebra ([3]{5j is provided by the Hilbert space H aux = 
L2[SU(2) 2n , dUf] where dUt is the Haar measure. The states have the form ip(Uf). The oper- 
ators Uf are diagonal and the operators Ef are the left invariant vector fields on each SU(2). The 
operators Et-i turn then out to be the right invariant vector fields. The operator associated to the 
volume Vt turns out to be the standard loop-quantum-gravity volume operator that is constructed 
in terms of Et. The states that solve the gauge constraint Q are labeled by SU(2) spin networks 
on the graph A*, which has a node for each tetrahedron and a link for each face of A n . A basis 
of these is given by states \ jf, it), where / = 1, 2n and t = 1, n range over the links and the 
nodes of the graph. These are defined by 

^ jfH (u f ) = (u f \j f ,Lt) = ®, nW(u f ) ■ ® t i t (io) 

where IL&(U) are the matrix elements of the spin-j representation of SU(2) and "•" indicates the 
contraction of the indices of these matrices with the indices of the intertwiners it dictated by the 
graph A*. For details, see j5]. 

With a scalar field, the Hilbert space becomes H aux = L2[SU(2) 2n , dUA (g) L2[R n ], with a 
(generalized) basis \jf, it, (f)t) and the states can be written in the form 

^{jf,H,<t>t) = {jf,kAt\^)- (11) 

In this basis the operator <p t is diagonal while p^ t = —i^~ t - 

If all constraints are first class, they can be quantized d la Dirac. The quantum Hamiltonian 
constraint can be defined in two alternative forms. The first, a la Thiemann, is obtained rewriting 
([7]) in the Thiemann's form 

C t = £ eV'f" Tr[U ff U^{U r ,V t }}~0 (12) 
ff'f& 



3 Ultralocal scalar fields are distinguished by the independent temporal development of the field at each spacial 
point. 
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and then denning the corresponding quantum operator by replacing the Poisson bracket with the 
commutator. Here the sum is over all ordered triples of distinct /'s and e" * is the parity of the 
ordered triple. The second possibility is to write directly the quantum operator which corresponds 
to the regularization of the Hamiltonian constraint used earlier in loop quantum gravity |22j 

C t = V t C t = Tr[Uff>E f Ef] ~ 0. (13) 
ff'& 

This form is more handy in the present context. Multiplying ([9| by Vt, we can rewrite the full 
quantum constraint equations in the form 

Sfi>= (~ P% + C t )i; = 0, (14) 

This set of n Hamiltonian constraints can be combined into a single one, introducing a "Lapse" 
N = {Nt} and writing 

S(N)ip = N t S t ip = 0, ViV. (15) 

t 

This concludes our definition of a family of finite-dimensional inhomogeneous quantum cosmolo- 
gies. We have one of these for each triangulation A n of a three-sphere. The hypothesis that we 
put forward is that they give an approximate description of the quantum behavior of our inhomo- 
geneous universe, increasingly accurate with n. This hypothesis can be seen as following naturally 
from the cosmological principle. 



IV. "DIPOLE" COSMOLOGY 

Consider the simple case obtained by taking n = 2 and the natural triangulation of the a three- 
sphere obtained by gluing two tetrahedra by all their faces. 4 A£ is then the graph formed by two 
nodes joined by four links 

A* 2 = 

The gravitational variables are (Uf,Ef),f = 1,2,3,4. We have two Hamiltonian constraints, one 
per each node, which we call C\ and C*2. Define 

C = V X C X + V 2 C 2 (16) 
and rewrite the constraints in the equivalent form 

S = C-+f(pi+4 2 )~0, (17) 

D = (VxC, - V 2 C 2 ) + -p)J ~ 0. (18) 

4 An interesting possibility, which we leave to the reader, is also to consider the n = 1 case defined by the graph 
Al — QO ' name ly a single tetrahedron with two couples of faces identified. 
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But using ([TJ and ([2]), it is easy to see that 

Vid = V 2 C 2 . 

Therefore D reads simply 

D = vl- PL ~ 0. 



(19) 



(20) 



which shows immediately that the Poisson bracket algebra between the two hamiltonian constraints 
closes. 

The gravitational Hilbert space is L2[SU(2) ] and a basis of spin network states that solve the 
gauge constraint is given by the states \jf, it) = |ji, j'2, jz, J4, ti, £2}- The action of one gravitational 
Hamiltonian constraint on a state gives 



C\jf,H) = ^ZC fr \j f ,L t ) 



(21) 



where, each term of the sum comes from one of the terms in the sum in / and /' in (13). More 
explicitly, we have 



e,5=±l 



(22) 



because the operator U12 = U\U 2 1 in (13) multiplies the terms n jl (f7i) and IP 2 (£/2) and 

un j (u) = u^ 2 (u)W(u) = c+ n j+1 / 2 (u) + c-W-^v). (23) 



The matrix elements Cj f ^ 2 can be computed with a straightforward exercise in recoupling theory 
from (13), and with some more algebra, from (12). In a different notation, in terms of the wave 
function components, we can write 



where Cj j f ^ t vanishes unless ej = for two and only two of the four j's. The scalar field variables 
are <fii, 4>2- Taking these into account leads to the wave functions ip(Jf, i n , fin), and (14) gives the 
dynamical equations 



k(3 2 

2 



2 + P)A,2 



J 2 

d 2 



v=o,±i 
d 2 



(25) 



(26) 



The coefficients C can be computed explicitly from recoupling theory. They vanish unless two ej's 
are zero. Equations fl25|26| ), defined on Hilbert space H 2 = L 2 [SU ~(2) 4 / 'SU \2) 2 \ <g> L 2 [R 2 } define a 
quantum cosmological model which is just one step out of homogeneity. 
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V. BORN-OPPENHEIMER APPROXIMATION AND LQC 



We now ask if and how LQC is contained in the model defined above. The state space H 2 
contains a subspace that could be identified as a homogeneous universe. This is the subspace 
H hom C Hi spanned by the states \j,j,j,j, ij, ij, 4>, <f>) where tj is the eigenstate of the volume that 
better approximates the volume of a classical tetrahedron whose triangles have area j. However, 
the dynamical equations 

the inhomogeneous degrees of freedom cannot remain sharply vanishing in quantum mechanics, 
due to Heisenberg uncertainty. Therefore it would be wrong to search for states that reproduce 
LQG exactly, within this model. In which sense then can a quantum homogeneous cosmology make 
sense? 

On the basis of the above discussion on the cosmological principle, the answer should be clear. 
The cosmological principle is the hypothesis that in the theory there is a regime where the inho- 
mogeneous degrees of freedom do not affect too much the dynamics of the homogeneous degrees 
of freedom, and that the state of the universe happens to be within such a regime. In other words, 
the homogeneous degrees of freedom can be treated as "heavy" degrees of freedom, in the sense of 
the Born-Oppenheimer approximation, and the inhomogeneous one can be treated as "light" ones. 
Let us therefore separate explicitly the two sets of degrees of freedom. This can be done as follows. 

First, change variables from the group variables Uf G SU(2) to algebra variables Af G su{2), 
defined by expAf = Uf. Following what is done in loop quantum cosmology [23J, let us fix a 
fiducial su{2) element ujf G su(2) for each face /. (This can be interpreted as the logarithm of the 
holonomy of the fiducial connection along the link dual to /.) We choose for simplicity a fiducial 
connection normalized as \uf\ = 1, and such that the four vectors u)f are normal to the faces of 
a regular tetrahedron centered at the origin of su(2) ~ i? 3 . Using this, we can decompose our 
variables as 



25|26 ) do not preserve this subspace. This is physically correct, because 



A f = cu f + a f , (27) 
E f = p uf + h f . (28) 

We need two conditions in order to fix this decomposition uniquely. First, we require that p is 
determined by the total volume 

V = ph (29) 

Second, we require that c is its conjugate variable, that is 

8ttG 

{c,p} = — = 1. (30) 

The variable c can then be identified with the corresponding variable using in quantum cosmology. 
We also define AV = V 2 -V 1 , so that V 1)2 = \{V± AV). 

Inserting the decomposition described above in the quantum Hamiltonian constraint ( 12 ) gives 



C t = - ^ [e a ^ +a fe~ cu f'~ a t'e~ cu t'~ a t' [e CUJ f" +a f" , V ± AV]] . (31) 

ff'fet 

Let us now decompose this constraint into two parts, the first of which depends only on the 
homogeneous variable c. This can be done keeping only the first term of the expansion of the 
exponentials in af and af, and only the V term in the volume term. That is, we write 
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C t = X -C hom + Cf (32) 

where 

c hom = £ r^j-gO^g-cw^g-^^^yjj = l^ftom (33) 
//7"6t 

The interpretation of this spilt is transparent: C hom gives the gravitational energy in the homoge- 
neous degree of freedom, while C\ n gives the sum of the energy in the inhomogeneous degrees of 
freedom and the interaction energy between the two sets of degrees of freedom. Finally, we write 
the homogeneous variable 4> = <pi + 4>2 and </>_ = <f>\ — 4>2- 

Following Born and Oppenheimer, let us now make the hypothesis that the state can be rewritten 
in the form 



rp(Uf, (j> t ) = iphom(c, cj)) ipi n h(c, 4>; a f , 4>- 



(34) 



where the variation of V'mft with respect to c and 4> can be neglected at first order. Here iphom 
represents the quantum state of the homogeneous cosmological variables, while tpi n h represents 
the quantum state of the inhomogeneous fluctuations over the homogeneous background (c,(j)). 
Inserting the Born-Oppenheimer ansatz (34) into (15), and taking N\ = N2, we have the equation 



Hi Hi (j^ ~ ~ ■ 

2 V'irah ^hom + ^ ^hom q ,2 ~~ V'inC ° m ' l Phom ~ C m iphom'4 ) in = 0. (35) 

Dividing by Vwt/W this gives 

I ggVW, _ C hom ^ hom = 2 W.^ inh + C inh jj hom ^nh (36) 

Whom Whom Winh WhomWinh 

Since the left hand side of this equation does not depend on the inhomogeneous variables, there 
must be a function p(c, 4>) such that 



k d 2 
2 



2 Whom - C hom Whom - PWhom = 0, (37) 



Hi d 2 C™^ 1 WhomWinh , /00 n 

o ^I^^inh H ; = PWinh- (38) 

2 d(j)i Whom 

The second equation is the Schrodinger equation for the inhomogeneous modes in the background 
homogeneous cosmology (c, </>), where p{c, <j)) plays the role of energy eigenvalue. The first equation 
is the quantum Friedmann equation for the homogeneous degrees of freedom (c, <j>), corrected by the 
energy density p(c,(p) of the inhomogeneous modes [21]. At the order zero of the approximation, 
where we disregard entirely the effect of the inhomogeneous modes on the homogeneous modes, we 
obtain 



1 °L 

2 d(t) 2 ' 



Whom = C^hom- (39) 



Let us now analyze the action of the operator C hom , defined in (33). Notice that c multiplies the 
generator of a U(l) subgroup of SU(2) 4 . Therefore it is a periodic variable c G [0,47r]. We can 
therefore expand the states Whom(c, </>) in Fourier sum 



11 



(40) 



The basis of states (c \ = e % ^ c l 2 satisfies 



p%) = (n/2)*\n) 
-4sin 2 (c/2)|^) = |/j + 2>-2|/i> + |/i-2) 



(41) 
(42) 



which we shall use below. The homogeneous Hamiltonian constraint (33) can be rewritten as 

c hom = £ Tr ^ cog |u + 2 sin |oy)(cos |l - 2 sin ^u^e^f" [e CUJ f" , V] 
ff'f" 



//'/" 



cos 2 1 1 + 2 sin | cos | (ay - ay/ ) + 4 sin 2 1 ay a// ) e~ CLJ /" [e^f" , V] . (43) 



Consider the action of the last factor on the state 



3 3 



Now observe that we can write 



3 9c/ \ 3 dc 

( / u/2)§e^ c/2 - e^f'k ((/i/2)l - ioy») 3 e ^/ 2 -n") 
((M/2)^l - (W2)l - e^ 2 . 

3 

((^/2)1 - iufr) 2 = a(/x)l + /3(/x)oy* 



(44) 



(45) 



where the coefficients a(/x) and /3(/x) can be easily computed squaring this equation. We write 
a(/i) = (/u/2) 2 — a(n). Bringing everything together 5 the only term that survives is 



<U//)+ ( 3^ 2§ j sin 2 | 



= [D + (fi) e +ic + £>-(^) e~ ic + £>V) ] eW2 

where D+(/z) = L>~(m) = - V^VsPif-i) and = ^a(fi) + 2\^v^/3(/i). 

Using (42), this gives 

C*H M > = Z>+(/z) |/x + 2) + Z>V) I/,) + D-( M ) | M -2). 



(46) 



The full equation ( 39 ) can be written as 
2 



2 *M) = C + (M) ^/i + 2» + C°( / u) *(/x^) + C-(/i) *(/x-2,, 



(47) 



(48) 



We are using Tr[tj/] = 0, Tr[ujfUjfi] = — ^e^e^qai, = — ~ cosOjfi and the symmetry of the homogenous case for 
which cos9ff> = cos8ffn = | and Tr [u)j uif =3*22 . 
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where we have written \I/ for tphom an d C =t '°(//) 



^Z>±.0( M ). 



kV2 



Equation (48) has the struc- 



ture of the LQC dynamical equation. Thus, LQC appears in the zero order Born-Oppenheimer 
approximation of a loop quantum gravity quantization of a finite number of degrees of freedom 
of the gravitational field, truncated according to the approximation dictated by the cosmological 
principle. 



VI. CONCLUSIONS AND PERSPECTIVES 

We leave a number of questions open. In particular: (i) the extension to real Immirzi parameter 
7; (ii) the inclusion of more realistic matter fields; (iii) understanding the relation between the 
ip{v,4>) homogeneous states and the full ip(jf,Lt,4>t) states in the spinnetwork basis; and (iv) the 
detailed comparison with the LQC quantization of the homogeneous universe, in particular in 
relation to the jl quantization scheme [1, 27 . The relation with more conventional cosmological 
perturbation expansion (see [LOj ) need also to be investigated. 

We observe that the discussion on the cosmological principle given in Section 2 and the observa- 
tion on the use of coarse triangulations in Section 3 bear also on the discussion on the viability of 
the recent computations of background independent re-point functions [19] . The conclusion of the 
discussion was the possibility of representing a large universe in terms of a coarse triangulation, 
and therefore, in the quantum theory, in terms of states based on graph with a small numbers of 
nodes. In cosmology, as in the calculation of the re-point functions, we observe that general rela- 
tivity admits an expansion in which a region of length scale L can be approximated by neglecting 
wavelengths much smaller than L, and described by spin networks with a small number of nodes in 
the quantum theory. In other words, the intuition that dynamics of nearly flat space can only be 
described using states with a high number of nodes is misleading. If this was the case, cosmology 
itself would be ill conceived. 

The idea of describing cosmological evolution using Regge calculus has been explored in the 
past. See for instance |l8j . Here we have adapted this idea to quantum cosmology, where it turns 
out to be particularly suitable for a loop quantization. The construction in [18] indicates that it is 
possible to have a 4d triangulation sliced by 3d triangulations equal to one another, and therefore 
suggests the possibility of writing a spinfoam version of the models introduced here. For instance, 
consider the kinematics of the A2 model and let ^4(jab> L a), 0, = 1, 5 be the vertex amplitude of a 
spinfoam model, as for instance that introduced in [23] . We can interpolate between two A2 with 
the triangulation A5 defined by the boundary of a four-simplex. In fact, collapsing four of the five 
tetrahedra of A5 into a single one gives precisely A2. This collapse is a 4-1 Pachner move, which 
can be realized interpolating a four-simplex. Therefore we can have a transition A2 — ► A2 via an 
intermediate A5. See Figure 1. 




FIG. 1: Spinfoam evolution of the dipole cosmology. 
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Accordingly, we can write a A2 — > A2 transition amplitude as the A2 — ► A5 — > A2 amplitude 
defined by the transition amplitude A(j a b, i a ) for each move. For instance 



-4-0'/, H]j'f, 4) = $jxj'i Yl A (3t,3c,j' t , l<c, tl, t'l)A(3t,jc,j't,jab, t>c, <-2, 4)- ( 49 ) 



where c = 1,2, 3. Repeating this step four times, over different points of the triangulation generates 
spacetime with the £3 x [0, 1] topology [18]. Is this dynamics related to the canonical one defined 
here? 

In summary, the models presented in Section 3 open a systematic way for describing the in- 
homogeneous degrees of freedom in quantum cosmology. In particular, they open the possibility 
of checking whether the bounce scenario that is characteristic of the homogeneous theory survives 
in a inhomogeneous context. The simplest possibility is to analyze the cosmological evolution in 



the simple A2 model described in Section 4, defined by the two equations (25) and (26). Do these 
equations govern semiclassical wave packets undergoing the cosmological bounce? If the answer 
to the above question is positive, the solution would also provide a concrete description of the 
fluctuating geometry at the bounce. The state ip(jf, it, fit = ^bounce) would give such a description 
explicitly. It is tempting to begin speculate on the possible cosmological role of the fluctuations of 
the inhomogeneous degrees of freedom at (or near) the bounce. Could they play a role in structure 
formation? For inflation? 

The derivation of the structure of the LQC dynamical equation presented in Section 5 sheds light 
on the relation between LQC and full loop quantum gravity. In particular, we have argued that 
the first should be searched as a Born-Oppenheimer approximation, as suggested by the hypothesis 
that goes under the name of cosmological principle. 

Finally, we point out the existence of the term p(c, <p) in the quantum Friedmann equation 



(28). Its physical interpretation is clear: it represents the back reaction of the quantum fluctuating 
inhomogeneous degrees of freedom on the dynamics of the scale factor. Again, it is tempting to 
begin to speculate on the possible cosmological role of this energy density. Does it play a role in 
structure formation? For inflation? In relation to the cosmological constant? 



We thank the audience of the ILQG seminar, where this work was first presented, for the useful 
questions and an enlightening discussion. We thank Bianca Dittrich, John Barrett and Simone 
Speziale for useful exchanges. 
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APPENDIX A: A SIMPLER HAMILTONIAN CONSTRAINT 



Instead of using the form (12) of the Hamiltonian constraint, we can also use the simpler 



densitized form (13) which corresponds to the regularization of the Hamiltonian constraint used 



earlier in loop quantum gravity in [22]. Here we show that this regularization gives a simpler 
dynamical equation, but with the same structure as in LQC. 



In the quantum theory, the decomposition (28) becomes 



9 r 

L f = u f - + L f . 



where LfC = 0. Inserting this decomposition into the Hamiltonian constraint (13) gives 



//'St 



dc 



(Al) 



(A2) 



As before, we us decompose this constraint in two parts, the first of which depends only on the 
homogeneous variable c. This can be done keeping only the first term of the expansion of the 
exponentials in aj and df, and only the term quadratic in in the derivative part. That is 



a 



This can be rewritten as 



horn 



Uf'Uf] 



ff'et 



d_d_ 

dc dc 



s lc 



horn 



(A3) 



Qhom = ^Tr[(cos|lL+2sin|c i ; / )(cos|lL-2sin|w / /)w / /LJ / ] 



//' 
1 



c 
2 L 



l dc 2 



o 2 



(-e tc + 1 - e~ lc )- 

The action of this operator on the states iphom(f J >i 4>) is therefore easily computed 

1, 



C h ° m ^horn(v,' 



(A4) 



(A5) 



Bringing everything together, the full equation (39) reads 



2 



C + (/i) IphamiV + 2, 4>) + C°(/i) IphomiV, 4>) + C (/i) IphomiV - 2, <f>) + Tr^^homiHi <f>) = 0. (A6) 



where the coefficient take the simple form C^(n) = — C°(/i) = We do not know if this simple 
equation gives the same phenomenology as the one used in LQC. 
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